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ABSTRACT 

In star formation, magnetic fields act as a cosmic angular momentum extractor that 
increases mass accretion rates onto protostars and in the process, creates spectacular 
outflows. However, recently it has been argued that this magnetic brake is so strong 
that early protostellar disks - the cradles of planet formation - cannot form. Our 
three-dimensional numerical simulations of the early stages of collapse (< 10 5 yr) of 
overdense star-forming clouds form early outflows and have magnetically regulated and 
rotationally dominated disks (inside 10 AU) with high accretion rates, despite the slip 
of the field through the mostly neutral gas. We find that in three dimensions, magnetic 
fields suppress gravitationally driven instabilities which would otherwise prevent young, 
well ordered disks from forming. Our simulations have surprising consequences for the 
early formation of disks, their density and temperature structure, the mechanism and 
structure of early outflows, the flash heating of dust grains through ambipolar diffusion, 
and the origin of planets and binary stars. 

Subject headings: binaries: general — circumstellar matter — methods: numerical - 
stars: formation — stars: general — stars: low-mass, brown dwarfs 



1. Introduction 

Over the past decade numerical simulations have enabled the exploration of the central physi- 
cal questions regarding the nature of star formation - from the collapse of an initial dense gaseous 
molecular cloud core to the formation of a star and its associated protostellar disk and outflow, 
and the emergence of a star's planetary system. The formation of disks and jets during star 



formation is central to many of these issues, but very little is known about the earliest phases 
of their evolution. How is the initial excessive angular momentum associa ted with the star's 
natal core removed - through magnetic braki ng (|Basu &: Mouschoviad 11994 ) and then o utflows 



(Ba nerjee k. Pudrita 120061 : iPudritz et al.1 120071 ). or by spiral density waves in disks (e.g. lLarson 
2009)? Do multiple stars form through gravit a tional fragmentation of cor es or massive disks 



( Hosking fc Whitworthl 12004 : IPrice fc Bate 



2007 



Hennebelle fc Tevssierll2008l )? What is the sig- 



nifica nce of outflows and jets that are launched b efore most of the mass has collapsed into the 
disk (jLynden-Belll 12003: Ba neriee k, Pudrita 120061) in a wide variety of young stellar systems - 
from brown dwarfs (IWhelan et al J 120051 : iMachida et al.ll2009i ) to massive stars (jArce et al.l 120071 : 
Baneriee & Pudrita 120071 )? 



Extraction of angular momentum by magneti c fields that thread the collapsing gas may be 
too efficient, according to recent two-dim ensional ([Mellon &: Lill2008l ) and three-dimensional ax- 
isymmetric (jHennebelle fc Fromangj l2008l ) simulations, preventing the formation of rotationally 
dominated disks, e ven when the effects of imperfect coupling of the field with the gas are included 
(jMellon &; Lill2009i ). These results seemingly contradict the observations which clearly show that 
disks are present around mos t if not all young stars, ev e n in environments in which th e magnetic 
field is expected to be strong (jHartmann &; CalvetJll995l ). iHennebelle &; Fromangj (|2008l ) performed 
three-dimensional ideal MHD simulations of collapsing cores using a barotropic equation of state, 
concluding t hat no rotati o nally dominant structure is formed from 10 to 100 AU for highly magne- 
tized cores. iMellon &: Lil (|2008l ) performed two-dimensional ideal MHD simulations on collapsing 
singular isothermal toroids using a barotropic equation of state and an inner boundary at 6.7 AU 
(effectively a sink particle), finding that even moderately magnetized disks could not form. Such 
two-dimensional models impose a high degree of mathematical symmetry and therefore miss the 
formation of bars and spiral waves in disks. 

We improve on previous results through this three-dimension al adaptive mesh re finement 
(AMR) investigation which includes a full treatment of the cooling ( Banerjee et alibood) and the 
finite coupling of the magnetic field to the pre-stellar gas (ambipolar diffusion, Puffin fc Pudrita 
20081 ). We find that ordered disk-like structures can emerge on scales < 10 AU at early times (< 10 5 
yr) in magnetized systems. We have omitted a sink particle in this study as they are expected to 
affect the solution to within a few sink radii. Without sinks simulations are indeed limited to early 
disks, although they offer a full solution to the region within 10 AU where the heated core forms. 



Methods 



We model an isolated star forming core wit h a rotating Bonnor-Ebert (BE) sphere ( Bonnor 



1956: lEbertlll955l ) commonly observed in nature (lAlves et alJl200ll ; lLada et al.ll2003l : iTeixeira et al 
20051 ). using the FLASH AMR code. The refinement criteria resolves the Jeans' length (Aj 



\Jnc 2 s /Gp) by at least eight cells, limited by the ambipolar diffusion time step. The ambipolar 
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diffusion time step is 



tad = 0.5- 



Ax 2 
VAD 



0.5 



Pol AD PiPn Ax 2 



(1) 



IAB 2 

where Ax is the cell length, t/ad is the diffusivity, 7ad = 3.28 x 10 1 3 g _1 cm 3 s _1 is the ion-neutral 
collisi onal coupling, B is the magnetic field strength, and pi = 2.32 x 10 -25 (p n /10~ 18 ) ' g cm -3 
(e.g. buffin fc Pudritzll2008l ). 



We construct an initial BE sphere with a radius of R = 6.9ro = 6. 9c Scorc / \J 4nGp c (where the 
critical radius is R c = 6.49ro) and add a 10% over-density on top with a 10%, m = 2 perturba- 
tion to break axisymmetry. Such perturbations model the buffeting that observed star-forming 
cores undergo from supersonic turbulent gas motions that prevail in the surrounding low density 
molecular gas (our perturbation accounts for a 0.04% change in gravitational binding energy, so 
this is a rather quiescent core). Our BE spheres have a mass of I.I8M0 and initial central densi- 
ties of p c = 6.07 x 10~ 18 g cm~ 3 . T he magnetic field is added to the core in the z direction - in 
agreement with recent observations (jKirk et al.ll2009l ) - such that the ratio of the thermal to the 
magnetic energy in the gas, known as the plasma beta (f3 = 2c 2 /v A = 46.01) is constant, where 
c s = 0.27 km s" 1 is the core's sound speed and the Alfvn speed is va = B/^A-Kp ~ 0.74c s /T (the 
latter relation for critical BE spheres). The mass-to-flux ratio T is a measure of the ratio of the 
gravitational to the magnetic energy in gravitating objects. If T < 1, the magnetic field is strong 
enough to prevent gravitational collapse and the core is called sub-critical, otherwise it is called 
supercritical. We set the mass-to-flux T = 2-kG 1 I 2 Y i / B = 3.5, where X is the colum n density, 
stem ming from observations of magnetic fields in early cores which indicate T ~ 1 — 4 (ICrutcher 



1999j) and confirmed through simulations of cloud-scale turbulence (jTillev fc Pudrita 120071 ). The 
rotation values of our m odels were taken from previous simulations of early star-forming clusters 
(jTillev <fe Pudrita 120071 ). We take the extremes and average values of rotation and interpret these 
as low, moderate and high rotating cores (Table [1]). The moderate rotation model set has a ratio 
of ro tational to gravitational binding energy of /5 ro t = 0.046, in line with observations (jLada et al 

20oah . 



For each model set, we run a non-magnetic (hydrodynamic), a perfectly coupled (ideal mag- 
netohydrodynamic) and more realistic, partially coupled (ambipolar diffusion) case for a total of 
nine simulations (summarized in Table [1]) . 



3. Do disks form? 

By the time the maximum surface densities in the moderate rotation model set reach a value 
of S c = 4.2 x 10 3 g cm -2 , 0.1 M & is contained within 156 AU (about 10% of the core's mass is 
contained within 0.003% of its initial volume). The mass of the protostar (taken to be all material 
inside 1 AU or 10~ 3 M©) is only 6 % of the total mass in the disk (taken to be everything inside 10 
AU or 10~ 2 M Q ). At this point, the non-magnetic, perfectly coupled and partially coupled clouds 
are 0.160, 0.194, and 0.175 million years old respectively (corresponding to 5.9, 7.2 and 6.5 free-fall 
times, where % = 0.027 Myr is estimated using p c ). An extended dense region has taken shape, 
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the gas has been heated, and in the cases where magnetism is present, disk winds have started. 



We present in Figure HJ snapshots of the structure of the region at the center of the collapse, 
taken at this time for our three different cases. Despite having identical initial conditions, the 
non-magnetized, moderate rotation case has formed two bars, driving off-axis rotational modes 
which themselves create an off-axis bar. The wobbling in the collapse occurs through accretion 
along shocks and density perturbations in the gas. This stands in stark contrast to the obvious 
disk-like structures and outflows present in both ambipolar diffusion (Figure 1(b) and ideal MHD 
cases (Figures 1(c) a nd (d)0. Strong oscillations of the disk radius were noted in the simulations of 



Price fc Bate! (|2007l ) which were also dominated by massive disks in their early phase. They also 



found that this behavior was suppressed by magnetic fields. 

Most importantly, the magnetic disks, by preventing strong bars from forming, are rotationally 
dominated whereas the uncoupled disk is not! This is illustrated in Figure Eta), where we see that 
the ambipolar diffusion case has a larger rotationally dominated disk (v^ > Wi n f a n) than the perfectly 
coupled case. We note also, by plotting f^/wmfaii of the most evolved perfectly coupled state of the 
collapse, that the rotational support is growing with time while hydrodynamic rotational support 
becomes increasingly unstable. In Figure Eta) we also plot surface averaged plots of the ratio of 
the rotational velocity to the Keplerian velocity - the rotational velocity of a gas in a stable orbit 
around an enclosed mass. We find peak ratios of about 0.1-0.5 in all cases of the moderate rotation 
model (in fact some pre-binned values exceed 1 in the ambipolar case). 

In the limit of low rotation, hydrodynamic spiral waves and bars in the purely hydrodynamic 
case are confined to a region of about 1 AU in radius. This has little effect on the braking of the 
cloud (Figure Etd) and rotationally dominated disks are allowed to form outside of 1 AU (Figure 
EJb). The hydrodynamic case is much more spherical in nature while the magnetized cases remain 
fairly flat, despite efficient magnetic braking. In the limit of low rotation, Figure 2(b) shows that 
the rotation in the hydrodynamic case is nearly Keplerian and the magnetic case is much more 
sub-Keplerian, whereas in the case of moderate rotation shown in Figure 2(a) the opposite is true. 
Thus in the extreme low rotation limit where spiral modes and bars are strongly suppressed, early 
hydrodynamic disks do form. 

Due to the bar, the distribution of specific angular momentum j{R) (where R is the cylindrical 
radius; Figure Etd) of the hydrodynamic case is comparable to the cases where magnetic braking is 
present. In the low rotation model set, the bar is small and j{R) does not compare as well to the 
magnetized cases. The hydrodynamic and perfectly coupled cases of the moderate rotation model 
set have similar maximal accretion rates of M accr ~ 4 x 10 -4 M yr" 1 , while the ambipolar case 
reaches maximal accretion rates of only M accr ~ lx 10 -4 M yr _1 due to weaker magnetic braking 
and a stable disk. In the low rotation model set, mass accretion rates are M accr ~ 2x 10 -4 M yr" 1 
for the magnetized cases and Maccr ~ 3 x 10~ 4 M yr _1 for the hydrodynamic case. 



1 We provide movies of Figures[TIa)-(c) online illustrating the density and magnetic field line structure on all scales 
in the simulation. 
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In the high rotation model set however, no cases form rotationally dominated disks as all three 
cases have formed some sort of bar on large scales (Figure [3]) . Clearly, it is through the suppression 
of instabilities in the collapse that magnetic fields promote early disk formation. Disk formation in 
systems without magnetic fields typically must wait until the mass of the central star dominates 
that of the surroundings sufficiently to suppress these instabilities. 



4. Early outflows and magnetic fields 

Figures QJb)-(d) are snapshots of the structure of outflows that are launched at these early 
times. A major result, shown by comparing Figures [TJb) and (c), is that outflows occur even when 
ambipolar diffusion is active. Outflows in a non-ideal MHD collapse have been demonstrated in two- 
dimensional ambipolar diffusion without drift heating ( Mello n &; Lill2009l ) and in three-dimensional 



using an ohmic diffusion approximation (jMachida et al 



2007). Our analytical understanding of the 



early outflow mechanism is the magnetic tower (|Lynden-Belll 120031 ) . In this picture, the toroidal 
field in the disk is generated through the rotating flow and confined there by the accretion ram 
pressure. As the field winds up, it reaches a critical strength in which the toroidal pressure is strong 
enough to overcome the ram pressure, moving wrapped field lines away from the disk and taking 
perfectly coupled gas with them - starting with the initial release in Figure QJc) and the more 
developed magnetic bubble shown in Figure QJd) . This acts to torque down the disk, removing 
angular momentum and increasing accretion rate! In the ^bipolar dnTusion picture/these fielj 
lines can seep through the ram pressure as they are not so tightly coupled to the accreting gas. 

We have also discovered the existence of a centralized, high-speed component to the ambipolar 
diffusion case's outflow (y > 0.4 km s -1 in an region roughly 1 AU in size). This component 
resembles a jet that is commonly observed in the evolved perfectly coupled collaps^f]. We note that 
the base of this component to the outflow lies on a cushion of heat generated by the friction between 
ions and neutrals as toroidal field lines seep through the accreting gas. This spike in drift heating 
i s typical in a mbipolar diffusion mediated C-shocks, regions where ion-neutral friction is strong 



(Wardlc 



1991 



This local heating effect produces strong gradients in temperature with T > 110 
K in regions roughly 5 AU in size, both above and below the midplane. Within these regions are 
features roughly 1 AU in size with T > 1000 K (not unheard of in C-shocks). By comparison, 
the perfectly coupled outflow has a sharp temperature shock at about 50 K leading its outflow. 
In contrast, the low rotation model set shows very little drift heating and a similar weak outflow 
(Figure G]Jc) . While the vertical field threading the disk B z is identical for all low and moderate 
rotation, magnetized cases (Figure [2(f), the toroidal field Bm in the low rotation, magnetized cases 
is nearly 2 orders of magnitude smaller than in the moderate rotation, magnetized cases. 

In Figure 0c) we plot the ratio of mass loss rates M w i n( j to mass accretion rates M accr for 



2 A movie of the development of Figure [T^d) from Figure QJc) is provided online. 

3 The velocity structure of the perfectly coupled outflow is seen in our online movie of the evolution of Figure [ljd) 
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the moderate and low rotation model sets. For comparison we also plot this ratio for the evolved 
outflow of the perfectly coupled case at later times (labeled 'evolved ideal'). The evolved outflow 
is very efficient in the moderate rotation, ideal MHD case - the central outflow component exceed s 
speeds of 5 km s _1 and observed outflow rates of M w ind/Afaccr ~ 0.1 (IHartmann fe Calvetill995l ). 
Unfortunately, we cannot fully develop the ambipolar outflow much beyond 10 5 yr due to time step 
constraints. However, from this early footprint, we are led to believe that the outflow mechanism 
is working. In the low rotation, ideal MHD case, the outflow rates are significantly diminished 
by nearly 2 orders of magnitude in comparison to mass accretion rates, less than what is typically 
observed. This m ay have consequen ces on feedback effects in star forming clusters and the assembly 
of massive stars (jWang et al.ll2009l ); protostellar energy injection is dependent on the rotational 
energy of the star forming core. 



5. Fragmentation and binaries 



In supercritical cores, three-dimensional ideal MHD calculations have shown that strong per- 
turbations (on the order of 50%) are needed in order for cores with r ~ 2 to frag ment and form 



binaries, even with high rotation (jPrice fc Bate 



2007 



Hennebelle fc Tevssierl 120081 ). We find that 



none of the low and moderate rotation cases, magnetized or not, fragment. We therefore ran a high 
rotation model set with four times the angul ar rotation than our ear lier model set (roughly 20 times 
the rotationa l energ y, an extreme rotation (jTillev Pudrita 120071 )). similar to the simulations of 
Price Bate! (J2007J) with a 10% perturbation amplitude. In Figure [3] we plot three-dimensional 
density contours for all three cases of the high rotation model at the largest common central column 
density (S c = 68 g cm -2 ). The hydrodynamic case quickly fragments into a binary with a wide 
separation of roughly 1000 AU. The perfectly magnetized case forms a bar and does not fragment, 
stabilized through its magnetic pressure as well as through the effects of magnetic braking. Mean- 
while, the ambipolar diffusion case produces an intermediate result. Early on a ring is formed which 
fragments into a bar-like object with a small companion at about 1000 AU. 

For sufficiently high rotation ambipolar diffusion allows binaries to form without the need of 
strong (e.g. > 10% in amplitude) perturbations. For moderate rotation rates however, fragmenta- 
tion can take place only if large perturbation amplitudes (~ 50%) are used in the initial conditions. 
Such a perturbation would amount to a 1% change in our BE sphere's grav itational binding energy. 
This relatively low level of 'turbulent' energy has been measured in B68 (|Lada et al.ll2003l ). This 
would suggest that a low level of turbulent energy (e.g. a few percent of the gravitational energy) 
is required in order to reproduce the overwhelming tendency for quiescent cores to form binaries in 
face of strong magnetic support. 
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6. Conclusions 



In the early stages of a purely hydrodynamic collapse of a moderately rotating core, material 
joins the protostar by accreting through a chaotic series of bars and spiral waves. Our results 
show that in magnetized collapses however, the magnetic field suppresses these wave modes, and a 
small, regular disk appears at the earliest times. The weakening of magnetic control by ambipolar 
diffusion is insufficient to guarantee the formation of binary stars in typical cores with moderate 
rotation. Our results suggest that modest turbulent amplitudes (> 10%) appear to be required. 
Regarding planet formation, the central density structure of early disks (Figure [2je) falls off as 
S oc r — (i- 7- 2 - 5 ) much more quickly with disk radius than do protoplanetary disk models at later 
times, wherein E oc r -1 ' 5 or r . As the collapse winds up the field early outflows appear - even 
for partially coupled disks - and feed angular momentum and mechanical energy back into the star 
forming neighborhood. The density and magnetic structure of ideal and partially coupled disks are 
quite similar - the main difference is in the strength of the wound up magnetic field (at least an 
order of magnitude weaker in the ambipolar diffusions case). 

Finally, we note that evidence for early localized intense drift heating near the disk at the 
accretion shock, which heats materials up to 1,000K and beyond in a localized region, may be 
preserved in the observed composition of comets and meteorites. The rapid heating, and subsequent 
cooling of those crystalline Mg-rich silicate materials that passed through the accretion shock and 
through the region of high ambipolar heating is reminiscent of heating events t hat must have 



occurred for some of these materials seen in cometary grains (IWooden et ajj 120071). These event s 
have been attributed to shock heating by spiral waves out to 10 AU in disks (jDesch et al.ll2005l ). 
As seen in Figure 1(b), only a portion of the accreted disk material will pass through this localized 
heated region. Subsequent radial turbulent mixing of this flash heated material with the bulk of 
material that passed through a more gentle heating environment could in principle contribute to 
the wide mixture of thermal histories preserved in cometary grain materials. 



7. Acknowledgements 

We thank Robi Banerjee, Eve C. Ostriker, and James Wadsley for fruitful discussions. D.D. was 
supported by McMaster University and R.E.P. by the Natural Sciences and Engineering Research 
Council of Canada. We are pleased to acknowledge the SHARCNET HPC Consortium for the use 
of its facilities at McMaster University. The software used in this work was in part developed by the 
DOE-supported ASC/Alliances Center for Astrophysical Thermonuclear Flashes at the University 
of Chicago. This research was supported in part by the National Science Foundation under grant 
PHY05-51164. 

REFERENCES 

Alves, J. F., Lada, C. J., & Lada, E. A. 2001, Nat, 409, 159 

7 



Arce, H. G., Shepherd, D., Gueth, F., Lee, C.-F., Bachiller, R., Rosen, A., k Beuther, H. 2007, in 
Protostars and Planets V, ed. B. Reipurth, D. Jewitt, k K. Keil, 245-260 

Banerjee, R., k Pudritz, R. E. 2006, ApJ, 641, 949 

-. 2007, ApJ, 660, 479 

Banerjee, R., Pudritz, R. E., k Anderson, D. W. 2006, MNRAS, 373, 1091 
Basu, S., k Mouschovias, T. C. 1994, ApJ, 432, 720 
Bonnor, W. B. 1956, MNRAS, 116, 351 
Crutcher, R. M. 1999, ApJ, 520, 706 

Desch, S. J., Ciesla, F. J., Hood, L. L., k Nakamoto, T. 2005, in Astronomical Society of the 
Pacific Conference Series, Vol. 341, Chondrites and the Protoplanetary Disk, ed. A. N. 
Krot, E. R. D. Scott, k B. Reipurth, 849-+ 

Duffin, D. F., k Pudritz, R. E. 2008, MNRAS, 391, 1659 

Ebert, R. 1955, Zeitschrift fur Astrophysik, 37, 217 

Hartmann, L., k Calvet, N. 1995, AJ, 109, 1846 

Hennebelle, P., k Fromang, S. 2008, A&A, 477, 9 

Hennebelle, P., k Teyssier, R. 2008, A&A, 477, 25 

Hosking, J. G., k Whitworth, A. P. 2004, MNRAS, 347, 1001 

Kirk, J. M., Crutcher, R. M., k Ward-Thompson, D. 2009, ApJ, 701, 1044 

Lada, C. J., Bergin, E. A., Alves, J. F., k Huard, T. L. 2003, ApJ, 586, 286 

Larson, R. B. 2009, ArXiv e-prints 

Lynden-Bell, D. 2003, MNRAS, 341, 1360 

Machida, M. N., Inutsuka, S.-i., k Matsumoto, T. 2007, ApJ, 670, 1198 
-. 2009, ApJ, 699, L157 

Mellon, R. R., k Li, Z.-Y. 2008, ApJ, 681, 1356 
-. 2009, ApJ, 698, 922 

Price, D. J., k Bate, M. R. 2007, MNRAS, 377, 77 

Pudritz, R. E., Ouyed, R., Fendt, C, k Brandenburg, A. 2007, in Protostars and Planets V, ed. 
B. Reipurth, D. Jewitt, k K. Keil, 277-294 

8 



Teixeira, P. S., Lada, C. J., & Alves, J. F. 2005, ApJ, 629, 276 
Tilley, D. A., & Pudritz, R. E. 2007, MNRAS, 382, 73 
Wang, P., Li, Z.-Y., Abel, T., k Nakamura, F. 2009, ArXiv e-prints 
Wardle, M. 1991, MNRAS, 251, 119 

Whelan, E. T., Ray, T. P., Bacciotti, F., Natta, A., Testi, L., & Randich, S. 2005, Nature, 435, 652 

Wooden, D., Desch, S., Harker, D., Gail, H.-P., & Keller, L. 2007, in Protostars and Planets V, ed. 
B. Reipurth, D. Jewitt, & K. Keil, 815-833 



This preprint was prepared with the AAS IATfjX macros v5.2. 



Model Set 


Low Rotation 


Moderate Rotation 


High Rotation 


Sltg = n v /3ir/32Gp 


0.1 


0.3 


1.2 


Angular rotation (fi, s -1 ) 


1.18 x 10~ 13 


3.52 x 10~ 13 


1.41 x 10~ 12 


Rotational beta (/3 ro t) 


0.0052 


0.046 


0.74 



Table 1: Model parameters. Each model set has a mass of M core = 1.18M , core temperature of 
T corc = 20 K (c Scorc = 0.267 km s" 1 ), external temperature of T cxt = 200 K (c Soxt = 0.845 km s" 1 ), 
T = 3.5, and (3 = 46.01 and is run in three cases: hydrodynamic, ideal MHD, and ambipolar 
diffusion. 
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Fig. 1. — Inner disk of the moderate rotation model set at central surface densities of £o = 
4.2 x 10 3 g cm -2 for (a) hydrodynamic, (b) ambipolar diffusion, (c) perfectly coupled cases, and (d) 

Black and gray surfaces 



2.58 x 10 10 cm" 3 ) 



the evolved state of the perfectly coupled case at S c = 8.72 x 10 7 g cm~ 2 . 
are p = 1CP 12 g cm~ 3 (n = 2.58 x 10 11 cm~ 3 ) and p = 1CP 13 g cm~ 3 (n 
respectively. Temperature contours at T = 85 K are shown in blue. Outflow velocities contours, 
shown in red, are 0.3 km s _1 in (c) and (d) and 0.01 — 0.4 km s _1 in b. Magnetic field lines are 
shown as yellow tubes. (Animations of this figure are available in the online journal.) 
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Fig. 2. — Azimuthally averaged plots (density weighted, see iBanerjee fc Pudrita I2006I ) vs. disk 
radius for ideal MHD, non-ideal MHD, and hydro cases. In all figures the model sets are shown 
at a common central surface density of S c = 4.2 x 10 3 g cm -2 (moderate or mod. rotation) and 
S c = 3.4 x 10 3 g cm -2 (low rotation), when viewed face on. Panels (a) and (b) show the toroidal 
velocity in units of infall velocity and Keplerian velocity for the moderate and low rotation model 
sets, respectively ("evolved" refers to hydro and ideal MHD end-states). The dotted horizontal 
lines represent a velocity ratio of 1. The dotted vertical lines represent a transition radius of 7 and 
6 AU, respectively. Panel (c) shows the outflow efficiency of the magnetized cases as compared to 
the more evolved ideal MHD case. Panel (d) shows specific angular momentum plotted against 
enclosed mass. The surface density distributions are shown in panel (e), while the magnetic field 
distributions are shown in panel (f). 
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Fig. 3. — Large scale fragmentation of the high rotation model set for (a) hydrodynamic, (b) 
ambipolar diffusion, and (c) ideal MHD cases at similar scales and orientations. Density contours 
are given in the legend, and magnetic field lines are drawn as yellow tubes. 
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